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ABSTRACT 

The Spacebome Laser Ranging System is a proposed short pulse laser on board an orbiting space- 
craft (1,2,3, 4). It measures the distances between the spacecraft and many laser retroreflectors 
(targets) deployed on the earth’s surface. The precision of these range measurements is assumed 
to be about ±2 cm (5). These measurements are then used together with the orbital dynamics 
of the spacecraft, to derive the intersite vector betwe en the laser ground targets. The errors associ- 
ated with this vector are on the order of 1 to 2 cm. The baseline distances to be determined range 
from 25 km to 1 200 km. By repeating the measurements of the intersite vector, strain and strain 
rate errors are estimated. These quantities are essential for crustal dynamic studies which include 
determination and monitoring of strain near seismic zones, land subsidence, and edifice building 
preceding volcanic eruptions. The realizable precision for intersite distance determination is esti- 
mated to be on the order of 0.5 cm at 300 km and about 1 .5 cm at 1 200 km. The corresponding 
inaccuracies for the intersite distances are iarger, that is 1 cm and 3.5 cm respectively. The corre- 
sponding precision in the vertical direction is 1 cm and 3 cm. The accuracies in the vertical direction 
which can be achieved are 3 cm and 10 cm. These values were obtained for a six day observing 
period with 50% cloud covwr. 

It is evident that such a system can also be used for geodetic surveys where such accuracies are gen- 
erally not needed. 
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PERFORMANCE ANALYSIS 
OF THE 

SPACEBORNE LASER RANGING SYSTEM 


I. INTRODUCTION 

Tlic result of computer simulations demonstrating the precision in determining intersite distances 
using a spaceborne Laser Ranging System aredescribed. Repeated determinations of intersite distance 
generate estimates of crustal strain and strain rate. The study of strain and the variation of strain with 
time are essential parameters for earthquake research (6, 7), It will be shown below that measurement 
periods of only a few days can yield very high precision measurements wliich, because of their accuracy 
and the speed with which they are obtained, can provide a new dimension to eartliquake study. 

II. SYSTEM DESCRIPTION 

Tin Spaceborne Laser Ranging System consists of an orbiting spacecraft carrying a pulsed laser 
distance measurement system that sequentially measures the distance to a number of retroreflector 
arrays on the ground. Figure 1 shows the general concept of the system. The spacecraft ranges to 
corner retroflectors on the Earth’s surface as it passes overhead. The proposed laser system consists 
of a Nd YAG Laser with a 200 picosecond pul.se length and a repetition rate of 10 pulses/sec. The 
RMS range uncertainty of a single pulse at 5 to 10 photo electrons is expected to be 1 to 2 cm with 
a bias of a few millimeters (8). The ground target (9) will consist of a small corner cube array of 
retroreflectors mounted on a pillar. 

As the first reflector of a ground network comes into view of the spacecraft, an acquisition 
procedure is initiated that is expected to take 10 to 15 seconds. The procedure consists of a search 
for the reflector based on a priori knowledge of the reflector’s location and the position of the 
spacecraft. After the acquisition of the first retroflector, the laser makes 20 to 30 range measure- 
ments in a 2 to 3 second period and then swings on to the next reflector, taking less than 0.5 sec- 
onds for this operation even for the most widely separated reflectors. The la-’^r dwells 2 to 3 
seconds on the second reflector making range measurements and then moves on to the next. No 
acquisition time is expected to be necessary for tlie second and subsequent reflectors because 








^ 9. ... 



the relative location of the spacecraft and ground network wilt be updated with the corrected a 
priori positions obtained during the acquisition of the first reflector. 

On any particular pass of the spacecraft ever the ground targets, the spaceborne laser will range 
in a preprogramed fashion to a given reflector approximately three times, each for a 2 to 3 second 
period; once at a low to medium elevation on approach, once at a high elevation, and once at a 
medium to low elevation on the way out. For a 1000 km altitude orbit, a pass of the satellite 
over the target area will last about 10 minutes whk h iinpiiv*^ *hat about 60 reflector arrays could 
be surveyed on every pass over the region. The lowest elevation at which measurements are to be 
made is about 20 degrees to minimize atmospheric refraction. Fig. 2 shows the sequence of events 
as the spaceborne laser passes over a network. 


III. INTERSITE DISTANCE ESTIMATION-ANALYSIS 

In what is to follow, all vectors are referenced to a common geocentric Earth fixed cartesian 
coordinate system. Furthermore, it is assumed that relatively small effects such as polar motion. 
Earth and ocean tides are properly modelled. S represents the vector position of a laser retroreflec- 
tor (laser target). Suppose that laser range measurements are obtained from m satellite passes, 
where Xj, i = 1, 2, ... ,m, is a six-dimensional epoch state for the i‘'' pass. The position of the 
satellite during the i^** pass at time r can be expressed as 

U.(r) = U(X.,r) (1) 

Where the function U is obtained by integrating the equations of motion with initial condi- 
tions provided by Xj from epoch time r. The measurement of the “round-trip” travel time of a 
pulse sent from the satellite to the laser target on the ground at time r when scaled by the speed of 
light is essentially a range measurement. Hence the fundamental measurement is considered to be 
the range (see Fig. 3), that is: 


rj(T) = U] (r) Uj(T) +STs-2S^Ui(r)] 
Where the observation rj is obtained during the i*** pass. 


( 2 ) 
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Figure 2. Sequence of Events as Spacecraft Passes Over a Network. 



Arrange ail observations into a column vector 


r = 


ri 

rj 


■••nJ 


(nXl) 


(3) 


From Equations (I) and (2) rcan be modelled by the non-linear equation 

r=F(Z,L) + c (4) 

The vector Z contains epoch values of the parameters to be estimated (e.g. satellite state, reflector 
positions, etc.), and vector L contains values of the unadjusted parameters which are assumed to 
be known constants in solving the regression equations (e.g. gravity coefficients, measurement 
biases, refraction errors, etc.), and e is the zero mean measurement noise vector. It is assumed that 
the elements of e are statistically independent. The final errors in the adjusted parameters can be 
decomposed into a component due to measurement noise, an alias component due to errors in the 
unadjusted parameters, (eg. the bias in the range measurement) and a component due to the errors 
in the a priori orbital state. 

The errors in the measurement represented by Equation 4 can be approximated by a first 
order Taylor series expansion about some current nominal parameter values of the vectors Z' and L'. 
Then to first order 


6f= f- F (Z' , L') = A (Z - Z') + B (L - L') + e 
= ASZ + B6L + e 


where 


AS 


9F(Z,L) 

. az 


z = z' 

L = L' 


(5) 


B = 

where Z' and L' are the nominal values. 


3F(Z,L) 


L 9L J 

Z = Z' 
L = L' 
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A and B arc the sensitivity matrices associated with the adjusted and unadjusted parameters. 
Given an unbiased a priori estimate 2o with error covariance iitutcix Wq, the weighted least 


squares iteration for estimating 2 is given by 


Z = Z' + 


A'^W A + Wi ■' A’^W"'«F+Wo"' (Zo -Z') 


where 


( 6 ) 


W s E (eeT ) 

On the first iteration. Z' is usually chosen to be Zq but on subsequent iterations, Z'= ^ 

A 

The estimated vector Z is subject to deviations in the a priori estimate, measurements, and 
unadjusted parameters, that is 

Z = Z + 6^ (7) 

where 

5^ - 6Zo 

A A 

6Zq : Error in Z due to deviation in a priori estimate 

A A 

52^: Error in Z due to deviation in measurement 

A A 

6Z^ : Error in Z due to deviation in the unadjusted parameters 

6Z„= [aTw-*a + w;'i-‘ w;' 6Z„ (8) 

5Z^= [A^W-'A + W;' )■' (AT;v-*e) 

6\= fA'^WA + W;']”' (A^W-'B^SL 

If it is assumed that and 6L are uncorrelated, then the overall covariance matrix 

associated with the estimated parameter vector ^ is 

E(5i6#'^) = ( A’^' W‘ A + W;* 1 

H-IA’^W-' A + W;']“‘ (A'^W-'B)Wl (BT" W-‘A)[A'^W-‘ A + W;* ]-‘ (9) 

where 

Wl =E(6L6LT) 
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In Equation 4 the parameter vector 2 is represented by 6ni epoch states and 3C laser target 


positions, that is 


X„ 


.Sg J (6m + 3fi) X I 


Of interest are the relative positions of the laser targets with respect to each other. A con- 


venient way of expressing this is in a tangent plane baseline coordinate svstem. This is obtained by 


differencing the station coordinates with respect to a master station 


D = QZ 


where D is the matrix consisting of all the iniersite vectors related to a master station (target) and Q 


is a transformation matrix. 


For example, for C laser targets with target number one being the master station and m satellite 


states 


O I I I -I I O 


O 10 -I 


01 10 10 


(3S -3) X (6m + 3fi) 


where 


I = Identity matrbe 
O = Null matrix 


This differencing is followed by a rotation so that the errors are expressed as an along base- 


line component (in the tangent plane), a cross baseline component, and a vertical component. 


D = RD = RC^ 










where 




1 

1 ± 1 


Rs 

0 

^ 1 

1 : 1 

^ 1 0 

■ 1 : 


0~ 

0" 1 



(3C-3) X(3fi-3) 


(13) 


For each baseline there is a different rotation inatnx Rg. 

Thus the covariance matrix in baseline coordinates is 

E(5D6DT ) s RQE (SZ 5p ) qT rT 
where E (8Z 8^^) is given by Eqn. (9). 

The covariance matrix E (6D6D^ ) gives a measure of the uncertainty of the intersite vector 
D. 

The repeated measurement of D gives the neasurement precision. 

If D, is the 1 St determination of the baseline vector over the survey period Tg and Dj is the 
second determination of the same baseline vector after a period T^ (the resurvey period) then 

(D, -Dj) = RQ(2, -Zj) (14) 

and 

8(D, -Dj) = RQ(5Z, -fiZj) (15) 

The aim here is to obtain the measure of precision on the baseline vector Dj and D 2 , that is 
E |5(Bi, -D2)5(D, -D2)T|. 

Using Equation (15) together with (7) and (8) yields this measure, that is: 

E |8(D, -D2)6(D, - D2)ij = E(N, N/) + E(N2 N 2 ’’) + (?, -P2)E(8L6 -P 2 )’' 


where 


E(N^N/) = RQ(A^Vr' A + Wq"* )^"‘ Q^rT 

P^ = RQfA^W* A + W(,-‘)f* (A’'\r*B),^ 

k = 1, 2 

Note that the precision of the baseline vector is in essence noise limited, since the sensitivity 
to the uncertainty in geopotential inherent in matrix P|^ is almost the same from one survey 


(16) 






period Tg to the next survey period. Only second order temporal effects (ie, d'ag, earth rotation, 
solar pressure, refraction...) still remain and influence the precision of the intersite vector determination. 

Tlie precision in the baseline component of (D,- D 2 ) is On . Tlie corresponding elongation error 
is given by 


= T* 


(17) 


where 

I Dl = magnitude of vector D on the baseline length. 

Baselines are frequently resurveyed a number of times. The estimated intersite distances can 
be fitted to a linear regression line ( 1 0) of the form 

d(t) = d„+cit (18) 

• 2 
The elongation rates d, determined by least squares, have a variance o^ which can be deduced 


from the least squares solution of Equation 18. That is 

,2 “ 


0.1 


n2 tf 
i=i‘ 




(19) 


where n is the number of measurements of the intersite distance of line d(t) and is the variance 
of d. If the measurements are made over time intervals Tg and repeated T,^ time intervals apart, 
then at the end of T (days, weeks, . . .) (See Figure 4) 


( 20 ) 
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where 


Tg = Survey period 

T^ = Resurvey, period 

T = Total measurement period, 
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Figure 4. Relationship Between Survey, Resurvcy and Total Experiment Periods. 

II 







From this, the elongation rate error for a given baseline can be deduced. From Equation 



IV. SIMULATION RESULTS 


A series of simulations have been performed of a survey of the States of California using a Space- 
borne Laser Ranging System. In the simulations approximately 1 50 laser targets are distributed over 
California at a separation of 50 km (Sec Fig 5). The simulations estimate tlie precision with whicli tlie 
baseline distance can be obtained in the presence of noise and bias of the laser system; perturbations 
of the spacecraft motion and errors in the refraction calculations. 

The orbit of the satellite is assumed to be circular at 1000 km and 50 degree inclination. A 
medium inclination orbit was chosen because it provides ground tracks across California in almost 
orthogonal directions (southwest to northeast and northwest to southeast) as shown in Fig. 6. 

This distribution of satellite ground tracks provides a strong geometric distribution of range measure- 
ments. In contrast, a polar orbit provides only north to south and south to north tracks and these 
provide strong geodetic ties in the north-south direction but only weak control in the east-west 
direction. 

The simulations have been conducted over survey intervals (T^) of 1 , 3, and 6 days respectively 
assuming 50% cloud cover that for the six day observation period reduces the number of success- 
fully observed tracks over the area from 36 to 1 8, For all the simulations the data on the observed 
tracks is assumed to be taken at the rate of 10 pulses per second with a noise of 2 cm and a bias of 
0.3 cm. The effect of errors in the gravity field on the motion of the satellite were accounted for 
by adopting the GEM-IO covariance model of the gravity field derived from satellite tracking and 
surface gravity data! II). The effects of solar radiation pressure and air drag on the satellite were 
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Figure 5. Distribution of Ground Targets for California Sim 



(DEGREES) 

e Ranging System 


assumed to be in error by a constant percentage in their estimation of their effect on the solutions. 
Based on operational experience, the radiation pressure error is assumed to be 30% of the nominal 
and the atmospheric drag 20% of the nominal. Tlie effect of atmospheric refraction errors were esti- 
mated through a two parameter (ie pressure and PTK gradient) model (12). In this model, the tem- 
perature and pressure are assumed known at a limited number of locations in the region and are used 
to develop an atmospheric model of the whole region from which the temperature, pressure and PTK 
gradients at each of the laser targets can be estimated. The pressure error of ±1.0 mbar chosen for the 
error analysis represents a realistic estimate for the pressure measurement accuracy at weather sta- 
tions. Walter Hoeline specified in (13) accuracy figures of ±0.2 mb for a mercury barometer and 
±0.5 mb for an aneroid barometer. Either of these instruments may be used at a weather station. In 
addition, an analysis of meteorological data from 48 weather stations in Southern California and 
Nevada (12) has shown that the pressure measurement residuals from a regression fit vary from 
±0.4 mb to ±0.9 mb. 

Figure 7 shows the baseline precision as u function of the baseline length. This measure of pre 
cision for a 50 km baseline is about 0.3 cm and increases to I cm for a 1200 km baseline. For base- 
line lengths up to 300 km, the precision is primarily dependent on system noise, but for longer base- 
lines, say from 400 km to 1200 km, the uncertainty in the geopotential becomes the dominant error 
source. For a system noise of 3 cm, the noise curve in Fig 7. will be shifted upward by a factor of 
1.5 for baselines less than 400 km and then the unadjusted parameters, predominantly the uncertainty 
in the Earth’s gravity field, dominate the precision for longer baselines. 

Figure 8 shows the improvement in the precision in the measurement of baseline distances 
which results for increasing the survey period Tg. An increase in the survey period from 3 to 6 days 
results in only a small improvement in the baseline precision. 

From the knowledge of the precision with wliich baseline determinations can be made, the elon- 
gation rate precision can be calculated from Equation 2 1 . 

Fig 9. shows the elongation rate precision for a 50 km baseline as a function of the measurement 
program period T. For the CiUculation of elongation rate precision, a survey period Tj of six days 
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Figure 7. Baseline Precision vs Baseline Distance 
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BASELINE PRECISION (cm) 
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Figure 8. Baseline Precision vs Baseline Length 
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ELONGATION RATE PRECISION (PER YEAR) 
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ASSUMPTIONS: 

• SURVEY PERIOD (T3): 6 DAYS 

• BASELINE DISTANCE (D); 50KM 

• BASELINE PRECISION; 0.25CM 




• Tp = 6 DAYS 
O Tp = 14 DAYS 
□ Tp = 30 DAYS 




Tp = RESURVEY PERIOD (DAYS) 

T = TIME PERIOD BETWEEN 1ST AND LAST MEASUREMENT ( YRS) 


2 4 6 8 

T (YRS) 

Figure 9. Elongation Rate Precision vs Total Measurement Periods. 
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was assumed and an interval between surveys, the resurvey period, , of 6, 14, and 30 days was 
chosen. As can be seen from this figure, the Spaceborne Laser Ranging System has the capability of 
determining elongation rates to a precision of better than 1 x 10“* strain per year over a 2’/2 to 3 
year period T. Furthermore, an order of magnitude improvement in the elongation rate precision can 
be achieved over a five year period by making continuous measurements over about a week every 
other month. The figure also shows that little improvement in the elongation rate determinations is 
achieved by increasing the measurement frequency (ie decreasing the resurvey period ). This 
means power for the laser system can be conserved and thus the system’s lifetime extended without 
degrading the measurement precision. Finally, Figure 9 shows that elongation rate measurements to 
better than 1 x 10“^ per year may be possible within a decade or less of measurements. 

V. CONCLUSION 

It is shown that the concept of a Spaceborne Laser Ranging System has the capability to (a) 
determine baselines to a precision of less than 1 cm over distances of up to 1000 km; and (b) deter- 
mine the elongation rate to a precision of better than I part in 10* per year during a period of 2 yrs. 
Such a system could provide a capability to observe the precursory geodetic motions believed to 
occur before large earthquakes. Indeed, established on a global scale, with survey areas around all 
major seismic zones, the Spaceborne Laser Ranging System could provide the first real probability 
for “capturing” a magnitude 7.5, and above, earthquake. 

In addition, general geodetic survey work can be performed accurately and very rapidly with 
minimal incumbrance by the terrain. 
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